Aging effects in the quantum dynamics of a dissipative free particle: non-ohmic case 
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We report new results related to the two-time dynamics of the coordinate of a quantum free 
particle, damped through its interaction with a fractal thermal bath (non-ohmic coupling ~ u s with 
0<<5<lorl<<5<2). When the particle is localized, its position does not age. When it undergoes 
anomalous diffusion, only its displacement may be defined. It is shown to be an aging variable. The 
finite temperature aging regime is self-similar. It is described by a scaling function of the ratio t w /r 
of the waiting time to the observation time, as characterized by an exponent directly linked to 8. 
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In complex out-of-equilibrium systems such as struc- 
tural glasses and spin glasses, the trend towards equi- 
librium is accompanied by aging effects: for instance 
the two-time correlation functions of certain out-of- 
equilibrium dynamic variables may not be invariant by 
time translation, even in the limit of a large waiting time 
(or age). The fluctuation-dissipation theorem, which is 
valid for dynamic variables at equilibrum, is not verified 
in this case. Actually the study of aging effects and of the 
related violation of the fluctuation-dissipation theorem 
(FDT) is a fundamental problem of the physics of dissi- 
pative out-of-equilibrium systems Q. Several aspects of 
this problem still stand as largely open questions, espe- 
cially in the case of quantum dissipative systems. 

In order to be able to discuss these questions at any 
temperature, one has to work within a quantum frame- 
work, the time scale Ti/kT playing a crucial role in the 
low-temperature dynamics. Since aging effects are also 
encountered, not only in complex systems, but also in 
simpler systems, nor disordered nor frustrated it is 
very natural, to begin with, to carry out the study of 
quantum aging in this latter type of systems. 

One archetype of a simple quantum dissipative system 
displaying aging effects is a particle coupled to a thermal 
bath but otherwise free, that is evolving in the absence of 
potential. Aging effects on the displacement correlation 
function and the corresponding violation of the quan- 
tum FDT have recently been discussed m for a specific 
model of dissipation, namely the so-called ohmic model. 
It corresponds to a particle undergoing standard quan- 
tum Brownian motion The damped ohmic equa- 
tion of motion can be given a non retarded form in the 
classical limit, in which diffusion is normal. However, 
the ohmic model does not allow to handle all the physi- 
cal situations of interest, for instance those in which the 
damped motion is described by a truly retarded equation 
even in the classical limit and in which either localization 
or anomalous diffusion phenomena are taking place. 



In this Letter, we examine aging effects in this sim- 
ple system using for the dissipation a versatile model 
able to generate various damped equations of motion, 
either instantaneous or retarded in the classical limit. 
The dissipation is introduced via a linear coupling of the 
particle to a thermal bath having a continuous distribu- 
tion of modes of bandwidth uj c (Caldeira and Leggett 
model (|, 0). The central ingredient of the model is 
the product of the density of modes of the bath g(oj) 
times the squared coupling constant |A(w)| 2 , a product 
assumed to vary as uj s at frequencies u> <C u> c . In the 
ohmic model, the dissipative exponent 8 is equal to 1. 
The algebraic cases < 6 < 1 and 6 > 1 are known, 
respectively, as the subohmic or superohmic dissipation 
models 

In the following we first show that for < 8 < 2 the 
particle velocity equilibrates at large times and does not 
age. Then, turning to the study of the particle coor- 
dinate, we discuss the domain of (8, T) parameters for 
which aging is taking place, or not. Actually the ques- 
tion is a non-trivial one. Indeed the two following prop- 
erties have been demonstrated |§-pr|. First, at T = for 

< S < 1, the particle is localized, in the sense that the 
mean square displacement Ax 2 (t) — ([x(t) — x(0)) 2 } tends 
towards a constant at infinite time. Second, at T = for 

1 < 8 < 2, and also at finite T for < 8 < 2, Ax 2 (t) 
diverges at infinite time, the diffusion being anomalous 
except at finite T for 8 = 1. As a result, we find clear-cut 
behaviors: as far as the coordinate is concerned, the lo- 
calized particle does not age, while the diffusing one ages. 
In this latter situation we provide analytic expressions for 
the effective temperature and the associated fluctuation- 
dissipation ratio allowing to write a modified FDT at 
finite temperature. We show in particular that the aging 
regime which emerges at large times is self-similar. 

In the Caldeira and Leggett model, the particle is cou- 
pled linearly to a set of harmonic oscillators in thermal 
equilibrium at temperature T. The Hamiltonian of the 
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particle-plus-bath system reads, in obvious notations, 
2m 



where A„ is a real coupling constant. For uj > 0, the 
quantity 2tt g(ui) \X(uj)\ 2 = \mhujK[uj) is modelized as 
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where 7 is a coupling frequency and f c a high-frequency 
cut-off function of typical width uj c . The definition of 
K(oj) is extended to uj < by imposing that it must be 
an even function of uj. The particle position operator 
obeys the retarded equation of motion 

r* 1 

x(t)+ / dt'Mt - t')x(t') = -x(ti)k(t - U) H F(t), 

Ju m 

(3) 

where ti denotes the initial time at which the coupling 
is switched on. In Eq. (||), the inverse Fourier trans- 
form k(t) of K(uj) plays the role of a memory kernel, and 
F(t) is a linear combination of bath operators, acting as 
a stationary random force of correlation function 



^$eK(uj)huj coth ^e"^ (4) 
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with $eK{uj) = \K(yj) and f3 = (keT)' 1 . 

It has been demonstrated that, for < 8 < 2 the total 
mass of the particle and of the bath oscillators diverges, 
while for 8 > 2 it remains finite and can be considered 
as a renormalized mass I [n|. For < 8 < 2 (<5 ^ 1), 
the average particle velocity {v(t)) relaxes towards zero 
at large t— U like (t — ti) 5 ~ 2 (the average is taken over the 
variables of the initial state of the particle and over the 
bath variables). For 8 = 1, the relaxation is exponential. 
The initial value (v(ti)) being forgotten for < 8 < 2, 
this situation may in this sense be qualified of ergodic. 
For 8 > 2, the dynamics is governed at large times by 
a kincmatical term involving the renormalized mass and 
(v(ti)) is never forgotten §Q. 

In the following, we limit ourselves to the ergodic case 
< 8 < 2. Then the velocity equilibrates at large times 
and does not age. This property, already obtained in 
the ohmic case || , thus generalizes to non-ohmic models 
with < 8 < 2. The two-time velocity correlation func- 
tion depends only on the time difference and can be com- 
puted via Fourier analysis and the Wiener-Khintchinc 
theorem: 
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With the modelization (||) for K(uj) and a lorentzian cut- 
off function f c = uj^/(uj 2 + uj 2 ), one has: 
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Then, setting 
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one may attempt to define the coordinate spectral den- 
sity as C xx (uj) = C vv (uj)/uj 2 . If convergent, the inte- 
gral (duj/2Tr) C xx (uj) represents (x 2 (t)), a quantity 
which must be independent of t. Checking the small-cj 
behavior of the integrand with the chosen modelization 
for the memory kernel, one sees that this is only possi- 
ble at T = for < 8 < 1. In this case, the particle 
is localized and it makes sense to define its position in 
an absolute way as x(t) = J_ v(t')dt' . The two-time 
position correlation function 



C xx (t,t>) = -({x(t),x(t')} + ), 



(8) 



where the symbol {., .}+ stands for the anticommutator, 
only depends on t — t' and does not age. In particular 
one has 
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that is, in terms of the integrated velocity correlation 
function D(t) = L duC vv (u), 



/oo 
D{t')dt' 
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The asymptotic behavior at large times of D(t) allows to 
characterize the relaxation of the mean square displace- 
ment towards Ax 2 (00). 

In other cases, that is at T = for I < 8 < 2 and at fi- 
nite T for < 8 < 2, the integral (duj/2n) C xx (uj) di- 
verges and (x 2 (t)) and C xx (t, t') are infinite. The particle 
diffuses. Since it is then no more possible to define an ab- 
solute position, we focus the interest on the displacement 
x(t) — x(to) (t > to). This quantity does not equilibrate 
with the bath, even at large times: it ages. The mean 
square displacement Ax 2 (t) diverges at infinite time and 
D(t) represents the time-dependent diffusion coefficient 
(in an extended sense when diffusion is anomalous, that 
is at T = for 1 < 8 < 2, and at finite T for < 8 < 1 
and I < 8 < 2). 

Before introducing the appropriate FDT violation fac- 
tor, let us study in details the behavior of D(t). Since 
< 8 < 2, one can restrict the study to the infinite bath 
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bandwidth limit uj c — ► oo in which case calculations are 
more simple. One has: 



D{t) = 7 
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At finite T, it is interesting to discuss on the same footing 
the classical counterpart of D(t), namely D cl (t) deduced 
from D(t) by replacing coth (/3hu>/2) by 2//3hu in for- 
mula (|ll|). Several important features of D(t) and D cl (t) 
can be obtained by contour integration. 

At T = 0, D(t) is found to be the sum of a pole con- 
tribution, which exists only for < 8 < 1, given by the 
oscillating function 



D[t) po i e 
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sin fit, 
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where f2 and A are known functions of 8 and 7, and of a 
cut contribution behaving at large times as a power-law, 

D(t) cut ~ — ( 7 <) 5 - 2 sin 3 ^r(2 - 6), (13) 
m7r 2 

where T denotes the Euler Gamma function. 

At finite T, D cl (t) is also found to be the sum of an 
oscillating function, which exists only for < 8 < 1, 



D cl (t) po ic 



with <f> = 7r<5/2(2 — 8) , and of a cut contribution behaving 
at large times as a power-law, 
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The behaviors of D(t) and D cl (t) at several different 
temperatures are illustrated on Fig. [I] for S = 0.5 and on 
Fig. H for 5 — 1.5. Interestingly enough, for any given 8, 
the curves corresponding to different bath temperatures 
do not intersect. Actually, it can be shown that, at any 
fixed time t, D(t), like D cl (t), is a monotonously increas- 
ing function of T. 

For times t <C tth (*th — h/2nkBT) and for any value 
of 8, the curves for D(t) at finite T nearly coincide with 
those at T — 0, as it should. 

At intermediate times, and for < S < 1, an oscillation 
due to the pole contribution takes place in D(t) (and 
also in D cl (t) but with a smaller amplitude). For certain 
values of 5 and T, this oscillation may even result in 
negative values of D(t) during a finite time interval. 

At large times t ^> 7 -1 ,ith and for any finite T, the 
curves of D(t) and D cl (t) join together: when < 5 < 1, 
D(t) describes a subdiffusive regime, and, when 6 > 1, 
a superdiffusive one. At T = 0, for < S < 1, D(t) 



describes the relaxation of Ax 2 (t) towards Ax 2 (00) and, 
for 1 < 8 < 2, a subdiffusive regime. 
1,2 
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FIG. 1. In full lines, D(t) plotted as a function of for 
8 = 0.5 and for bath temperatures T = 0, fcsT = Tiy^n, 
ksT = ^7/77 (at T = 0, D(t) is not a diffusion coefficient, 
but characterizes the relaxation of Ax 2 (t) towards its limit 
value). In dashed lines, the corresponding D cl (t). 
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FIG. 2. In full lines, the quantum diffusion coefficient 
D(t) plotted as a function of yt for 8 — 1.5 and for bath 
temperatures T = 0, ksT — hy/2-K, ksT = hy/n. In dashed 
lines, the corresponding D cl (t). 

These features of D(t), especially the large times ones, 
are essential for describing the aging properties of the 
displacement correlation function, as defined by 

C xx (t, t';t ) = ~({[x(t) - aj(*o)], W) - x(t )}}+). (16) 
In the classical modified FDT can be written 
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as §| 

Xxx (t, = /?9(t - t')* cl ft, t'; t ) dC ^f' *°) , (17) 

where x^fAO i s the displacement response func- 
tion. For a diffusing particle, the fluctuation-dissipation 
ratio X c \t,t';t ) can be obtained from D c1 (t) and 
D cl (t w ) ||, where r = t — t' denotes the observation time 
and t w = t' — to the waiting time: 



X c \r,t w ) = 



D c \t) 



D c1 {t) + D*\t w )' 



(18) 



For any r and t w , one can define an effective inverse tem- 
perature as /?eff( r ) t w ) = /3A cl (r, Since X cl does not 
depend on T, the bath temperature is rescaled by a fac- 
tor 1/X cl larger than 1, due to those fluctuations of the 
particle displacement which take place during the wait- 
ing time. At large times (r,t w ^> 7~ 1 ,£th)> on e can use 
in formula (|l^) the asymptotic expressions of D c1 (t) and 
D c \t w ) as given by Eq. (|l|). Eq. (|lS|) then displays 
the fact that, in a subohmic or superohmic model of ex- 
ponent (5(0<(5<lorl<(5<2),a self-similar aging 
regime takes place at large times, as pictured by 



or (possibly anomalously) diffusing ones in which it only 
makes sense to consider the displacement, which displays 
aging (T = 0, 1 < S < 2, or T finite, < 6 < 2). The ag- 
ing regime at finite T is properly described by the large 
times expression of the classical fluctuation-dissipation 
ratio, which for 5 1 is a self-similar function of t w /r, 
solely parametrized by 5. 
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Interestingly enough, X cl ' ag and T 



•cl,ag 
eff 



are functions of t w /r, solely parametrized by 6. They do 
not depend on the other parameters of the model (i.e. 7 
or oj c ). For (5 = 1, one retrieves the results X cl ' ag = 1/2 



and T 



■ci,ag 



off 
cl,a 
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and T c g are algebraic functions of t w /r 



For any other value of i5, X cl,ag 
The limits 
t — > 00 and t w — > 00 do not commute. 

In the quantum case, the effective temperature T e g = 
(/c/3 e ff) _1 can be obtained from the following equation ||: 



Dt*(t) = D(r)+D(t w ). 



(20) 



Eq. @ also allows to define T eS at T = for 1 < 6 < 2. 
Since D(t) is a monotonously increasing function of T, 
Eq. ( p0| ) yields for T c g (r, t w ) a uniquely defined value. 

The curves representing (3 c s(r,t w ) as a function of r 
for S — 0.5 and S = 1.5 at a given finite temperature and 
for a given t w 3> 7 -1 , ith are plotted on Fig. ^. Quantum 
effects do not persist beyond times r ~ f tll . Thus, for 
times r > 7 ~Vth, X cl < a s(r, i„) as given by formula Jl9| ) 
allows for a proper description of aging effects. 

In conclusion, we have shown that the two-time dy- 
namics of the coordinate of a quantum dissipative free 
particle coupled to a thermal bath (coupling ~ uj s with 
< S < 2) displays extremely rich behaviors. According 
to the values of T and S, one may find either a local- 
ized regime in which the position can be defined in an 
absolute way and does not age (T = 0, < 5 < 1), 



FIG. 3. In full lines, the effective inverse temperature /3 e ft, 
as computed from Eq. (|2o|), plotted as a function of for a 
bath temperature ksT = ?ry /2n and two different values of S. 
In dashed lines, the corresponding classical effective inverse 
temperature /3^ ag , as deduced from Eq. @. 
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